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We explore the A/ = i rule and e'/e in K —>^ tttt tran- 
sitions using a Dyson-Schwinger equation model. Exploiting 
the feature that QCD penguin operators direct Kg transitions 
through 0"*"^ intermediate states, we find an explanation of 
the enhancement of A' — > 7r7r/=o transitions in the contribu- 
tion of a light cr- meson. This mechanism also affects e' /e. 



I. INTRODUCTION 

The A/ = 2 rule is an empirical observation: the 
widths for nonleptonic decays of kaons and hyperons that 
change isospin by one-half unit are significantly larger 
than those for other K and A transitions; e.g. 



^K°^{TnT)/^K+ 



660. 



In terms of the amplitudes A/^o^^+^- and Mj^o^^o^o 
that describe Kg — > vrvr transitions, the pure isospin-zero 
and isospin-two n n final states are described by 

Ao = ^(2M^o^,+,- + Mko^^o^o) , (2) 

A2 = -^{Mko^^+^- - M^o^^o^o) , (3) 

and the ratio in Eq. (n|) corresponds to: 

l/w := Re(Ao)/Re(A2) « 22 . (4) 

The analogous amplitude ratio for S'-wave A -^ t:N tran- 
sitions is 1^1/2/^3/2! ~ 80. 

The processes involved are nonleptonic weak decays so 
one necessarily encounters QCD effects in their analysis 
and the operator product expansion (OPE) can therefore 
be employed to good effect. Using the OPE the ampli- 
tude for a given transition is expressed as the expecta- 
tion value of an effective Hamiltonian: A — (Ticff) = 
X]i ai(/j,) (Qiilt)) , where /z is a renormalisation point. 
The coefficients: ai{fj,), are calculable in perturbation 
theory and describe short-distance effects. However, 
the expectation values of the local effective operators: 
{Qi{lj)), contain the effects of bound state structure; i.e., 
long-distance QCD effects, and must be calculated non- 
perturbatively. 

The transitions of interest herein are mediated by non- 
leptonic strangeness changing (AS* = 1) effective opera- 
tors. The simplest: 



?2 = s^Oj^Ui UjO'^dj, 



(5) 
(6) 



with O^ = 7^(1±75) and colour indices: i,j = 1, . 
have the flavour structure of the standard weak four- 
fermion current-current interaction, and there are eight 
other terms representing the QCD and electroweak (ew) 
penguin operators.R At least some of these must have 
large expectation values if the A/ = i rule is to be un- 
derstood. 

Another quantity that may be much influenced by the 
AS* = 1 effective interaction is the ratio e'/e. The indirect 
CP violating parameter: 



A{Kl -> mri^o)/AiKs 



TTTT/^O) 



(7) 



measures the admixture of CP-even state in K^ '■ for e 



(1) 0, CPli^. 



L/S) 



T\K 



L/S) 



i.e., they are CP eigenstates. 



e appears to be primarily determined by short-distance 
contributions from the weak nonleptonic AS" = 2 effec- 
tive interaction [g[. In contrast, e' measures the phase of 
the heavy-quark CKM matrix elements in the standard 
model and 
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(8) 



with |e| = 0.002280, experimentally.^ A nonzero value 
of e'/e entails direct transitions between CP-even and 
CP-odd eigenstates. e' is sensitive to the same penguin 
operators that contribute to the A/ = ^ rule, and hence 
is likely to receive significant long-distance contributions. 
The current generation of experiments |^ appears to be 
consistent and an average value of the ratio is B 



Re ( - ) = (2.1 ±0.46) X 10" 



(9) 



Qi results from QCD corrections to the weak current- 
current vertex. The penguin operators are generated by QCD 
and ew corrections but their flavour structure is different. 

^In Eq. (pi) we follow contemporary practice and make ex- 
plicit the TTTT-scattering phase shifts: (5o,2, in factoring out the 
phase $£/ = ^ -\-52 — So- Then, using the experimental obser- 
vations; So « 37°, 62 « -7°, ^e ^ f , one has ^^, - <E>e ^ 
and the imaginary part in Eq. (H) relates only to an explicit 
CP violating phase. 



A standard form of the AS* = 1 effective interaction at 
a renornialisation scale /i = 1 GeV is 

10 

GFj2c,{fi)Qdp), (10) 



n_/AS=l 



where Gp = Gf K*s Kd/V2, Ci{fi) = z^{^j,) + Tyi{n), 
T = -{VtlVtd)/{V:,Vud), and Kd,. • ■ , are the CKM 
matrix elements. (Direct CP violation is a measure of 
Im(T).) The coefficients: Ci{fi), at next-to-leading order 
are quoted in Ref. |g|, as are the operators: Qj. W e repro- 
duce the coefficients in the appendix, Eq. (All), but not 
the operators and note only that Q3a,5,6 are the QCD 
penguin operators; e.g.. 



q—u,d.s 

and Q7,8,9,io £Li"C the ew penguin operators; e.g.. 



Q8 = 



- s O^d 



q—u,d.s 



GnQjO+qi 



(12) 



where Cq is the quark's electric charge (in units of the 
positron charge) . The expectation value of the operators 
in Eq. (|lO|); i.e., the long-distance contributions, are the 
primary source of theoretical uncertainty in the estima- 
tion oiw and e'/e, Eqs. (|) and §). 

Herein we calculate the expectation values of the op- 
erators in Eq. ( |lO| ) using the Dyson-Schwinger equation 
(DSE) model of Ref. 0. The DSE framework is reviewed 
in Ref. |^ with some of the phenomenological applica- 
tions described in Refs. |7|-jri[|. It treats mesons as bound 
states of a dressed-quark and -antiquark with Bethe- 
Salpeter amplitudes describing their internal structure, 
and has already been used to explore CP violation in 
hadrons |]12|. We describe the calculation and its ele- 
ments in Sect. ||, and present and discuss our results in 
Sect. ID, Section |^ is a brief recapitulation. 



II. OPERATOR EXPECTATION VALUES 

A. Charged kaon decay 

The impulse approximation to the meson-meson tran- 
sitions mediated by Ti^g'^^ is straightforward to evalu- 
ate; e.g., in the absence of ew penguins only the operators 
Qi,2 contribute to K^ -^ 7r+7r° transitions and 

(7r+(pi)^°(p2)|gi|if+(p)) = ;^ E KUpi,P2), (13) 
T^{puP2)^iV2irZj 0-x.{k2\~hP2.-\p2) (14) 

X 2trZ2 / O^XK{ki;p2,Pi)^Aku-Pi)Su{ki) , 

iT2{pi,P2)^2V2tYZi / o;;xAk2;~lP2,-^P2) 

Jki Jk2 

X Oj^ XK{ki;p2,Pi)^Tv{ki;-pi) Su{ki) , (15) 



with the trace over Dirac indices only, and 

xAk; ^1, ^2) = Su{k + li) r,(fc; £1 + £2) Su{k - £2) , (16) 
XKik;£i,£2) = Ss{k + £i)TK{k;£i + £2) Su{k - £2) . (17) 

Here we use a Euclidean formulation with {7^,7,^} = 

2'^M''' 7^ =7m andp-g = X;LiP»9i- /fe '■= I d'^k/{2Tr)* 
is a mnemonic representing a translationally invariant 
regularisation of the integral, with A the regularisation 
mass-scale that is removed (A -^ cxd) as the final stage of 
any calculation, and Z2{fi, A) is the quark wave function 
renornialisation constant. Sf^u,s are the dressed-quark 
propagators (we assume isospin symmetry) and TH=K,Tr 
are the meson Bethe-Salpeter amplitudes, both of which 



(11) we discuss in detail in Sect. II B 



Using the Fierz rearrangement property: 

tr[0;GiO;G2] = -tr [0;Gi] tr[0;G2] , (18) 

where Gi_2 are any Dirac matrices, it is clear that Ti ex 
T2. Furthermore, the analysis for Q2 is similar and the 
result is the same so that 



(^+(pi)7r°(p2)|(ciQl + C2Q2)\K+{p)) 
C\ +C2 



V2 



N,{N, + l)Ti{pi,p2). 



This can be simplified using Q 

.A 
/ttPm = - \/2 iVctrZa / O^ X7r{k;-lp,- \p) 

Jk 

and 

-{P + Pi),f+\pl)~P2,f''"{pl) = 
2NctrZ2 iO^XK{ki;p2,Pi)^Ak 

Jki 



(19) 



(20) 



Pi)Su{ki), (21) 



where f^ are the Kg-^ semileptonic transition form fac- 
tors, so that 



{^+{p,)i.'{P2m%'='\K+{p)) 



Nr 



V2N, 

Ml{pi,P2) 



Gf (Cl +C2)Mlipi,P2); 



u P2.{p+pi)f+ (pD+plf'^'ipl) 



fw{m 



frhmi) 



(22) 

(23) 
(24) 

-1.0 



where the last line follows from ||1^ 

and mlf^ (— "t,^) ~ 0. 

We can compare our result with the contemporary phe- 
nomenological approach to K — > tttt decays, which em- 
ploys a parametrisation of Al 1 : 



Ml 



fTvirriK 



mi)B\ 



(3/2) 



(25) 



with the parameter B[ ' ' fixed by fitting the experi- 
mental width. One historical means of estimating J\4i 
is to employ the vacuum saturation Ansatz, which yields 

^(3/2) ^ ^ j^ .g ^j^^g ^jg^^ j^^^ g^g ^j^ ^^^ ^<j 

that our impulse approximation is equivalent to this 
Ansatz.^ However, agreement with the experimental 



value of Tj^ 



requires Bj^ 



(3/2) 



as can be seen 



using Eq. ( |A12| ). Thus, while the impulse approxima- 
tion is reliable for estimating the order-of-magnitude, it 
appears that an accurate result requires additional con- 
tributionsJj However, our primary goal is to identify a 
plausible mechanism for an enhancement of tttt/^q tran- 
sitions and this level of accuracy is sufficient for that pur- 
pose. Hence we proceed by adopting the contemporary 
artifice and use 



Mlipi,P2) 



U{mji 



mi)B\ 



(3/2) 



r(3/2) _ 1 
^1 - 2 

j(3/2) 



(26) 



In doing this we bypass the calculation of _Bj^ , which 
our elucidation of the impulse approximation has identi- 
fied as a real challenge for models whose basis is kindred 
to ours, and also for other approaches. 



B. Propagators and Bethe-Salpeter amplitudes 

Although the matrix element discussed above was ex- 
pressed in terms of dressed u- and s-quark propagators, 
and TT- and ii'-meson Bethe-Salpeter amplitudes, we ob- 
tained a model independent result without introducing 
specific forms. That is an helpful but uncommon simplifi- 
cation only encountered before in the study of anomalous 
processes such as tt" -^ 77 |1^ . 

In general, as reviewed in Ref. g, these quantities can 
be obtained as solutions of the quark DSE and meson 
Bethe-Salpeter equation. However, the successful study 
of an extensive range of low- and high-energy light- and 
heavy-quark phenomena ||,0-|l2[ has led to the devel- 
opment of efhcacious algebraic parametrisations. These 
were employed in Ref. |q] and we also use them herein. 

The dressed-quark propagator is 



CT^(x) = 2 m/ T{2{x + mj)) 



(27) 
(28) 
(29) 



^This is an exact algebraic constraint, which has been over- 
looked by other authors and consequently violated in fitting 
Tj^+^^+^o; e.g., Ref. Ilq] and references therein. 

*The impulse approximation, with dressed-propagators and 
vertices, has proven quantitatively reliable in many other ap- 
plications; e.g., Refs. [p|,|7|-|ll{, and in some cases corrections 
have been calculated and shown to be small |14,P_5| . The new 
feature herein is that the calculation is not self-contained; i.e., 
we presently rely on external input: the d in Eq. (lid). 



+J'{b(x) T{b{x) \bl + b{T{ex) 



aUx) 



with J^{y) 



1 



x + m 



f 



[l-T{2{x + rh}))] , 



(30) 



(1 - e y)/y, X = _pV-^^ "V = iTT-f/\ 
ag{x) = X(7g{p^) (7y{x) = \^ (Jy{p'^). Thc mass-scale, 
A = 0.566 GeV, and parameter valuesjj 



m 



bo 61 



u 0.00948 0.131 2.94 0.733 0.185 
s 0.210 0.105 3.18 0.858 0.185 



(31) 



were fixed in a least-squares fit to light- and heavy-meson 
observables 0], with these dimensionless u, s current- 
quark masses corresponding to 



to},^^^ = 5.4 MeV, mi^^^^ = 119MeV. 



(32) 



This algebraic parametrisation combines the effects of 
confinement and dynamical chiral symmetry breaking 
with free-particle behaviour at large spacelike p^ Q . 

The dominant component of the tt- and ii'-meson 
Bethe-Salpeter amplitudes is primarily determined by 
the axial- vector Ward-Takahashi identity ||l^,Q : 



where Bh ■= B^\"!:'^l and g 

fej =0.204, 6^ = 0.319; 



(33) 



(34) 



i.e., Bh is the quark quark mass function obtained from 
Eqs. (p^)-(pn) with ifif = and 6q replaced by the values 
indicated. With these dressed-propagators and Bethe- 
Salpeter amplitudes one obtains (in GeV) 



and (qq) 



Calc. 
Obs. [1] 

IGcV 



u 



Ik 



rriK 



0.146 0.130 0.178 0.449 
0.131 0.138 0.160 0.496 



(35) 



(0.220 GeV)3 



C. Neutral kaon Decay 

We now turn to thc transitions K^ -^ 
In comparison with K'^ — > tt+tt" there is a significant 
qualitative difference: all effective operators contribute 
to these transitions and furthermore the QCD penguin 
operators: Qs.e, and ew penguin operators: Qt^Sj can 



^e — 10~* in Eq. (EOl) acts only to decouple the large- and 
intermediate-p^ domains. The study used Landau gauge be- 
cause it is a fixed point of the QCD renormalisation group 
and Z2 ~ 1, even nonperturbatively [[13|. 



direct the transition through 0"'""'" intermediate states. 
This can be important because; e.g., Ref. |]1^ reports 
evidence of a broad scalar resonance in r — > z^T-Tr^Tr^vr" 
decays: 



mo++ ~ 1.12 mx , r, 



0+ + - 



0.54 GeV, 



(36) 



and with mQ++ w mn such a resonance could provide 
a significant contribution to the nonleptonic K'^ decays. 
We explore this hypothesis by allowing such a contribu- 
tion in our analysis: Kg —^ CTo++ —>■ tttt. 

Before proceeding further we note that a light scalar 
meson is a feature of DSE studies using a well-constrained 
rainbow-ladder truncation |[l^-pl[|. However, in the 0++ 
channel this truncation is very likely unreliable P2| , a fact 
we expect is entangled with the phenomenological diffi- 
culties encountered in understanding the composition of 
scalar resonances below 1.4 GeV |23]. While the lack of a 
reliable DSE truncation for scalar mesons prevents an ac- 
curate calculation of their Bethe-Salpeter amplitude and 
mass, they are nevertheless describable by such ampli- 
tudes, which herein we parametrise as 



T^{k;p) ^ Id^y 



1 



{kycoD 



2\2 



(37) 



where Id — 74 and tUg- is a width parameter. T„ is nor- 
malised canonically and consistent with the impulse ap- 
proximation {q± = q ± ■^p) 



Pp. = Nc tr 



A r 



r^iq;~p)^i^T4q;p)Siq-) 



T^{q]-p)S{q+)ra{q;p) 



dpi, 
dS{q- 



dpij 



(38) 



p-=-"i„ 



We separate the Qg contribution to the Kg -^ tttt tran- 
sition into two parts and consider first the new class of 
contributions, which introduce the a intermediate state: 

{7r{p,)n{p2)\Qe\K°ip)) = (39) 

{7r{p,}n{p2)\a{p)) D,{p') {a{p)\Q,\K"{p)) , 

where we represent a propagation by 

D,ip^) = l/[p^ + ml] (40) 

and employ the impulse approximation for the cttttt cou- 
pling 

M^^^{pi,P2) := {tt{pi)tt{p2)\(j{p)) = (41) 

2iV,tr / T,{k;p)Su{k++) 

xir^(fco+; -pi) 5„(fc+_) ir^(fc_o; -P2) Suik — ) , 

ka/B = k+{a/2)pi + {P/2)p2, which provides the basis for 
the calculation of ro++^7r7r- This combination of simple- 
pole propagator plus impulse approximation coupling to 



the dominant decay channel is phenomenologically effi- 
cacious; e.g., Refs. [|lO| . 
In impulse approximation 



{a{p)\Qe\K''{p))^V2N^x 



tr Zi 



iXK{ki;\p,\p) Ot, X'y{k2; 



h, 



(42) 



with Z4(/i,A) the mass renormalisation constant and 
Xa(^;^ij^2) an obvious analogue of x-n{k',iiT(-2) in 
Eq. (|l|). Using 

tr[GiO+G20;] - 2tr[Gi(l - 75)] tr[G2(l + 75)] (43) 

this yields 

-^^{a{p)\Q,\K\p)) = 

iV2NctTZj ii5XKiki;b,lp)\ 



V2NctTZ. 



4 / x^(fc2; 

k2 



ipr 



■h) 



(44) 



From Refs. [g^ 24| we identify the first parenthesised term 
as the residue of the kaon pole in the pseudoscalar vertex: 



irK := 



V2 Nc tiZi / 75 XK {ki ; \p, \p) 

Jki 



-f 2 



K 



(45) 



The second term is the scalar meson analogue in 
the scalar vertex but the vector Ward-Takahashi iden- 
tity, which is relevant in this case, does not make 
possible an algebraic simplification. The integral 
and its /i-dependence must therefore be calculated. 
That is straightforward when the renormalisation-group- 
improvcd rainbow-ladder truncation is accurate; e.g., 
Refs. |liyi3| , |l9|l , but not yet for scalar mesons. This is 
where the simple Ansatz of Eq. ( p7|) is useful: it yields a 
finite integral and we therefore suppress Z^ to obtain 



Ua{p)\Q,\K^{p)) 



V2 



(46) 



.A 

r_R-V2iVctr / Xa{k2\-\p,-\p) =■ TK ra{p^) ■ 

Jk, 



The result for Qc, is similar, but suppressed by a factor 
of 1/Nc, and the contribution of the ew penguins: Qy.s, 
can be obtained similarly. 

The other class of contributions, which do not involve 
a 0+"'' intermediate state, can be evaluated following the 
explicit example of Qi presented above. Only two addi- 
tional three-point functions arise: 



G^{P1,P2) 
-GkAP1:P2) 



(7r(pi)Tr{p2)\{uu + dd)\0) , (47) 

{7r-ip,)\su\K''{p)). (48) 



They are the scalar pion form factor and the scalar Ktt 
transition form factor, respectively, and Qkt;{pi,P2) can 



be expressed without additional calculation in terms of 
the Kp-i form factors [fuj: 



QkApi,P2) 



(49) 
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Pi 



9 

P2 



rud 



pI) 



ni~P2) + 



pI 



pI 



Pi 



f^i-pi) 



a result which follows from the vector Ward-Takahashi 
identity. A preliminary result is available [E5[ for 
Q^{pi,P2), which takes the form anticipated from cur- 
rent algebra. That is to be expected because correctly 
truncated DSE models provide a good description of chi- 
ral symmetry and its dynamical breakdown, as illustrated 
in a study of nn scattering ]2q| . This makes a calculation 
of Q^{pi,p2) unnecessary for our present study because 
we can adopt the form ^ [(rf ) = 3.76 GeV-^]: 

g'APi,P2) - -4^ [1 „ 1 (,^^)2 (p^ +p,)2] . (50) 

Jtt 

The matrix elements for the K — > tttt transitions can 
all be written 



M, 



K- 



^qcd 



'^cm^y^K^TTW ' 



(51) 



with the explicit forms given in the appendix and the 
pure isospin amplitudes defined in Eqs. (|), (^). 



III. RESULTS AND DISCUSSION 

Everything required for our calculation of the widths is 
now specified. There are two parameters: lUo- in Eq. (37); 
and nia in Eq (|4(]|). We determine them in a least-squares 
fit to: Tj^o^^+^-, r^o^^o^o, taken from Ref. ||l[; and 



rcr^(7r7r) in Eq. (36), and obtain (in GeV) 





Obs. 


Calc. 


mc 


1.12mK 


l.lAmK 


UJcr 




0.611 


^cr^(Tnr) 


0.54 


0.54 




5.055 ±0.025 


5.16 


-pXlO ^■' 


2.305 ±0.023 


2.11 


-pXlO"^-' 


0.0112 ±0.0001 


0.0116 



(52) 



which is a relative error on fitted quantities of < 4%.H 
This value of uja- corresponds to an intrinsic cr-meson 
I/wct, which is 0.84 r^; i.e, 84% of that of the 



size: 



J — 



p-meson determined in Ref. |5[| 



"We used Gf = 1.166 x 10"" GeV"", Vts = 0.0385, Vtd = 
0.0085, Vus = 0.220, 14d = 0.975, ImiVt'sYtd) = 0.000133, and 
a obtained from Eq. (All|). With r„{k;p) oc exp(-fc^/a;^) 
instead of Eq. (|^), m^ = 1.12 rriK, (^^ = 0.694 GeV yields 
exactly the same results for the calculated quantities. 



The widths in Eq. ( p2| ) are obtained from the calcu- 
lated amplitudes (in GeV with niK from Eq. ([3^)) 

\Mko^^+^- I = 2.7 X lO"'^ = 5.9 X lO"'^ TUK , (53) 

IM^o^^o^o I = 2.4 X 10"^ = 5.4 X lO"'^ tuk , (54) 

\Mk+^^+^o I = 1.8 X 10"** = 4.0 X 10"^ ruK ■ (55) 

For the pure isospin amplitudes we find (in GeV): 

Re(ylo) = 31.7 x 10"^ , Rc(yl2) = 1.47 x 10"** , (56) 

which are consistent with recent lattice estimates p8| and 
yield 



l/w = 21.6. 



(57) 



Our analysis also yields values of the parameters: 
bI ' ''^ ' , used in phenomenological analyses to express 

the operator expectation values 0. Of course, B\ — 
0.5, as discussed in connection with Eq. (26) and, using 



the formulae in the appendix, we obtain algebraically: 



(1/2) 



(1/2) 



BX'^> = B. 

We also calculate 



= B. 



(1/2) 



= B 



(1/2) 



= B. 



'.^m 



= 51 



(3/2) 



B, 



(1/2) 



B, 



(1/2) 



1.43 +(17.9), 



(58) 



where the second term is the contribution of the cr-meson. 
The non-CT contribution is large because of the strength 
of the ATtt transition form factor. If the vacuum satura- 
tion Ansatz is used to estimate the operator expectation 
values they are all = 1. That method does not admit a 
cr-meson contribution. 

Eliminating the ew penguin contributions yields a 
< 1% reduction in 1/w, which is consistent with the 
the magnitude of cxcm- Suppressing instead the cr-meson 
contributi on, while not affecting '^ k+^tt+tt'^ of course 
(see Eq. (^), yields r^o_^+^- = 1.3 x IQ-i^GeV, 
r^fo^^o^o = 1.1 X 10"" GeV, and \/w = 2.9. 

The value of e'/e follows from Eq. (H). Suppressing 
the cr-meson and ew penguin contributions we obtain 
e'/e = 128 x 10~^, which is ~ 60-times larger than the 
experimental average in Eq. (^). Including the cr-meson 
we find 31.3 x 10~^. To understand these results we note 
that Eq. (H) can be written 



w ImA 



V2 |e| RcAn 



1- 



1 luiAo 



IrciA 



(59) 



which makes clear that the ratio 
Im(Ao)/Re(A o) un less Im(A2) ^ 0. 



is determined by 
Noting that ci^2 



are real, Eq. (All), then it follows from Eq. (Al) that 
Ini(j42) = in the absence of ew effects. Hence our calcu- 
lated results are large because the pre-factor in Eq. ( [59| ) 
is large. The dependence on the a contribution is easily 
understood. The pre-factor is oc Im(Ao)/Re(Ao)^, which 
is large in the absence of the a contribution even though 
Im(Ao) and Re(Ao) are individually small. The a contri- 
bution adds simultaneously to Im(Ao) and Re(Ao) with 



a magnitude ^ 100-times larger than the original val- 
ues. Hence the final ratio is sensitive only to the relative 
strength of the a contributions, which is determined by 
the coefficients cs.g- 

Including both the a and ew penguin contributions we 
obtain 



e7e = 31.7x 10^ 



(60) 



from which it is clear that the ew penguins are a cor- 
rection of order aem as one would naively expect. In 
this case Im(A2) y^ 0. However, as observed above, the 
CT-meson enhancement responsible for the A/ — ^ rule af- 
fects the real and imaginary parts of Ag simultaneously 
so that (l/w) ImA2/ImAo remains negligible. 

If we employ the artifice of an ad hoc suppression of 
the a contribution to Im(Ao) while retaining it in Re(Ao); 
i.e., make the replacement 



c, Ms -^ Re(c, M3), i^ 5,6,7,1 



in Eqs. (^)-(^, we find 

e'/e = 2.7x 10"^ 



(61) 



(62) 



those recently inferred [|I8[, and the analysis is not sen- 
sitive to details of the model Bethe-Salpeter amplitude. 
However, this same mechanism yields a value of e'/e that 
is ~ 15-times larger than the average of contemporary 
experimental results unless a means is found to suppress 
its contribution to 1tii{Aq). 

If a light scalar resonance exists it will contribute in the 
manner we have elucidated and should be incorporated 
in any treatment oi K ^ tttt. 
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This artifice is implicit in the phenomenological analyses 
reviewed in Ref. ||j and that is why Eq. ( p2^ reproduces 
their order-of-magnitude. The small value is only possi- 
ble because in this case Im(Ao) is not cr-enhanced and is 
therefore of the same magnitude as Imyl2 /w ex acm/w, 
due to the 1/w enhancement factor. That factor survives 
because Re(Ao) is still magnified as required in order to 
satisfy the A/ = -^ rule. Currently we cannot justify this 
procedure.n 



IV. EPILOGUE 

We have demonstrated that estimating the K — > 7r7r/=2 
matrix element using the impulse approximation is al- 
gebraically equivalent to using the vacuum saturation 
Ansatz and yields a result that is ~ 2-times too large. 
The identification of a compensating mechanism that can 
correct for this overestimate is a contemporary challenge. 

We have also shown that the contribution of a light 
scalar meson mediated by the QCD penguin operators: 
Q^fi, is a plausible candidate for the long-range mecha- 
nism underlying the enhancement of if — > tttt/^o transi- 
tions .n A good description of that enhancement requires 
a mass and width for this 0++ resonance that agree with 



APPENDIX: COLLECTED FORMULAE 

The matrix elements for the K — j tttt transitions are 
all of the form in Eq. ( pl| ) with 



^'^K+^.+.o = ;^ Gj. (1 + ^) (ci + C2) A^i , 

"§-..+ .- - Gf {[c2 + C4 + ^(Ci + C3)]Xl 
+ 2(^C5-HC6)(A^2 + ^A^3)} 



M 



(Al) 
(A2) 



MSfn'^ 



ifO^TrOTrO 



C4 - Cl 



f (C2-C3)]>ll 



Gf[[ 

+2 (^C5 + ce) (M2 + 75-^3)} , (A3) 






T2^^{ 



Cl + ±cs (A4) 

-(l + ^)(c9 + cio)]A^i+3(^c7-fc8)A^^} , 

Gf [{^^CQ + cw)Mi (A5) 

(^C7 + cs) {M^2 + 2X2 + 71 ^fa)} , 






-^R-O-^^o^o — G} 



{[c7 + i;C8-(l+ 2^)C9 



(A6) 



i(l + ^Jcio] Ml - (^C7 + cs){M2 + 73X3)} , 



''NB. If this procedure is followed then mu 7^ nid isospin 
symmetry breaking effects also contribute significantly to e'/e. 
Q5,6 mediated scalar diquark transitions: (us) ;^ — > 
(Md)Q+°, are the s — > t-channel interchange of the interac- 
tion that herein produces the a-meson. They are a viable 



candidate for the mechanism that produces the AJ = i en- 
hancement for baryons. This was explored in Ref. [29], how- 
ever, the requirement therein that diquarks also explain the 
enhancement for mesons appears unnecessarily cumbersome. 



A^2 

M2 



rKQ^{pi,P2). Ml^r^gf,^ipi,P2) 



Ml 



M 



2 ' 



M3 = rKrcr (p-) D„ {p ) Mg.^^ (pi , P2) 



(A7) 
(A8) 
(A9) 



83, 22 (1999); P. Debu, sem inar at CERN, 18/June/99, 



withp^ = {pi+P2y 



-m^.pl 



pI 



These for- 



mulae make clear the operators that would be suppressed 
if Nc were large. Note that 



M' 



qcd 



o-M^To'.--M],'o'o). (AlO) 



This is not true of the complete amplitude. 

In our calculations we use values of the coefficients that 
correspond to our choice of Aqcd ~ 0.2 GeV: c,; = z,; + 

TVi, T 



-{v,:vu)/{v:,Vud)with§ 





Zi 


Vi 


1 


-0.407 


0.0 


2 


1.204 


0.0 


3 


0.007 


0.023 


4 


-0.022 


-0.046 


5 


0.006 


0.004 


6 


-0.022 


-0.076 


7 


0.003 


-0.033 


8 


0.008 


0.121 


9 


0.007 


-1.479 


10 


-0.005 


0.540 



(All) 



Using the alternative set listed in Ref. H then, with 
mo- = 1.06 mif and Wo- = 0.670 GeV, we obtain re- 
sults that differ from those in Eq. (152) by < 1%, and 
e'/e = 69.0 x 10~^ primarily because 2/5 in the alterna- 
tive set is 2.6-times as large. 

Fro m t h e co mplete matrix elements: Eqs. (pL) and 
Eqs. (Al)-(A6), we obtain the widths: 






K°- 



^K'^^ttOtt" —C{mK)\M 1(0 ^^o^o\ , 



C{x) = 



1 



1 



4ml 



167r x 
while the matrix element of Eq. ( [4l| ) features in 

r^^^^^-, = ^C(mcr) \MoTT7r{ml;Tnl,ml)\'^ . 



^(tttt) 



(A12) 
(A13) 
(A14) 

(A15) 



(A16) 
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